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In second-order scalar-tensor theories we study how the Vainshtein mechanism works in a spheri- 
cally symmetric background with a matter source. In the presence of the field couphng F{<f)) = e"^'''^ 
with the Ricci scalar R we generally derive the Vainshtein radius within which the General Rela- 
tivistic behavior is recovered even for the couphng Q of the order of unity. Our analysis covers the 
^SJ ■ models such as the extended Galileon and Brans-Dicke theories with a dilatonic field self-interaction. 

T-H ' We show that, if these models are responsible for the cosmic acceleration today, the corrections to 

, gravitational potentials are generally small enough to be compatible with local gravity constraints. 

(N- 

^ . I. INTRODUCTION 

• The observational discovery of the late-time cosmic acceleration [l| poses one of the most serious problems in modern 
PsJ I cosmology. Within the framework of General Relativity (GR) it is possible to realize the accelerated expansion of 
the Universe by taking into account a "dark" component of a matter source. The representative model of this class is 
^ ■ quintessence, in which the source of dark energy comes from the potential energy of a scalar field (/) [2]. However, it is 
Qi generally difficult to accommodate an extremely tiny mass required for the cosmic acceleration today (m^ 10^^'^ eV) 
^ • in the framework of particle physics f§\ . 

{^jrjl An alternative approach to the dark energy problem is the modification of gravity at large distances . One of the 
I— simplest examples is the so-called f{R) gravity, in which the Lagrangian / is a general function of the Ricci scalar 
rvi '• ^ M; functions f{R) are required to be carefully designed to satisfy cosmological and local gravity constraints 

^ ■ |6h1C|. In viable dark energy models based on f{R) gravity the mass of a scalar degree of freedom is large in high- 

■ density regions, so that the chameleon mechanism can be at work to suppress the propagation of the fifth force. 
^\ . Nevertheless, in the cosmological context, such models are plagued by the fine tuning of initial conditions associated 

■ with the oscillating mode of field perturbations 0, M, Tl3| . This property generally persists for the chameleon 
l models of dark energy which are designed to pass both cosmological and local gravity constraints [fl [l5| . 

■ Another representative model of dark energy based on the modification of gravity is the Dvali-Gabadadze-Porrati 
^— I ] (DGP) braneworld scenario in which the cosmic acceleration is realized by a gravitational leakage to the extra 
^-H . dimension. Although this model contains a ghost mode [17] in addition to the incompatibility with observational 

" data [18,] , it has a nice feature to recover GR in a local region through the Vainshtein mechanism [19]. This property 
J> . comes from the field self- interaction of the form (V0)^n^, which appears as a mixture of the transverse graviton 
■ '~j ■ with a brane-bending mode . The non-linear field interaction suppresses the propagation of the fifth force for the 
rN '. distance smaller than the so-called Vainshtein radius ry. 

^ ■ The Vainshtein mechanism was originally proposed in the context of a Lorentz-invariant massive spin-2 Pauli-Fierz 
■ " " " theory [2lj. The quadratic Pauli-Fierz theory possesses the van Dam-Veltman-Zakharov (vDVZ) discontinuity [12] 
with which the linearized GR is not recovered in the limit that the mass of the graviton is zero. Vainshtein showed 
that in the nonlinear version of the Pauli-Fierz theory there is a well-behaved expansion valid within a radius ry [3 ■ 
Although the nonlinearities that cure the vDVZ discontinuity problem typically give rise to the so-called Boulware- 
Deser ghost in massive gravity [2^, it is possible to construct nonlinear massive gravitational theories free from the 
ghost problem in the decoupling limit [24l]. Recently the Vainshtein mechanism was applied to the (new) massive 
gravity models [2^ and also to Galileon models [26, 27,]. 

In Galileon gravity the field Lagrangian is constructed to satisfy the Galilean symmetry d^cp ^ d^cf) + bf^ in the 
limit of flat spacetime [2^. The nonlinear field self-interaction XO(f>, where X = — (V0)^/2, appears as one of those 
terms [29J. The cosmology based on Galileon gravity has been extensively studied recently in the context of dark 
energy [2^ [s^, [3l| and inflation [s^] ■ For the covariant Galileon there exists a de Sitter solution with a constant field 
velocity. Thus the cosmic acceleration can be driven by the field kinetic energy without a potential. Moreover there 
are some viable parameter spaces in which the ghosts and Laplacian instabilities are absent [soj . 

Galileon gravity can be viewed as one of the specific theories having second-order field equations. The general action 
of scalar-tensor theories with second-order equations was first derived by Horndeski ^] in the context of Lovelock 
gravity. This issue was recently revisited by Deffayet et al. [34] as an extension of Galileon gravity (see also Ref. [35|). 
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The general _D-dimensional action derived in Ref. fsS] reproduces the Horndeski's action in four dimensions [36']. The 
Gahleon term XDcf), for example, can be promoted to the form G(0, X)n0, where G is an arbitrary function with 
respect to (f) and X. In fact, the dynamics of dark energy in the presence of the term G{(l), X)0(j) have been studied 
by a number of authors [ST^ . 

The Horndeski's action involves the Lagrangians of the forms £2 ~ Pi't': ^3 = ^G{(j>, X)D(p, £4 — 6*4(0, X)R + 
Gi.x X (field derivatives), and £5 ~ G5{<j),X)G^i,^{\i"-'\i"'(j)) + (Gs^x/S) x (field derivatives), where R is the Ricci 
scalar. G^i/ is the Einstein tensor, and Gi^x = dGi/dX (i = 4,5). For the functions G4 that depend on alone, 
e.g., G4 — F(0)/2, the Lagrangian £4 reduces to £4 = F{(f))R/2. For the choice F = Mpj, where Mpi is the reduced 
Planck mass, £4 corresponds to the Einstein-Hilbert Lagrangian. 

In this paper we shall study the Vainshtein mechanism in a spherically symmetric background for the Horndeski's 
second-order theories with G4 = F{(f))/2 and G5 = 0. We do not take into account the effects of the term G5 as well 
as the X-dependence in G4, but the presence of the nonlinear field interaction G(0, X)[2(j) can allow us to understand 
how the Vainshtein mechanism works in the presence of the nonminimal coupling F[(f>)R/2. Our analysis covers a 
wide range of gravitational theories such as (extended) Galileon, dilaton gravity, and Brans-Dicke theories with the 
nonlinear field interaction. 

This paper is organized as follows. In Sec. |TT] we derive the equations of motion for the action ([T]) given below in 
a spherically symmetric background. Under certain approximations we also obtain the simplified equations for the 
field as well as the gravitational potentials. In Sec. Illll we clarify how the Vainshtein mechanism works in general and 
constrain the forms of the action. In Sec. IIVI the general results derived in Sec. IIIII are applied to specific models. 
We also estimate the corrections to the gravitational potentials coming from the modification of gravity. Sec. |V] is 
devoted to conclusions. 



II. FIELD EQUATIONS IN A SPHERICALLY SYMMETRIC BACKGROUND 

We start with the following action 



S 



F{(f)) R + P{(t>, X) - G(0, x)n(t> 



+ Smig^i^: *m) : 



(1) 



where g is the determinant of the metric g^^, R is the Ricci scalar, F((/)) is a function of the scalar field 0, P{(j),X) 
and G{(j),X) are functions of and X — —g^^^d^(f)dv(j)/2, and Sm is the matter action. We assume that the matter 
fields ^'„j do not have direct couplings with the field 0. 

We derive the equations of motion in the spherically symmetric background with the line element 



(2) 



where ^(r) and <I>(r) are functions with respect to the distance r from the center of symmetry. For the matter action 
Sm we consider a perfect fiuid with the energy-momentum tensor — diag {—pm, Pm, Pm, Pm)- The (00), (11), (22) 
components of the equations of motion derived from the action ([l} are given, respectively, by 
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(3) 
(4) 



+ F" + — - <t>'^G,^ - 2(0" - $'0')^G,x = e'*(Pm + P) , (5) 
r 



where X — —e 2*0'2/2, a prime represents the derivative with respect to r, and a comma corresponds to the partial 
derivative in terms of or X (e.g., G,0 — dG/dcj)). The equation of motion for the field is 



2e-"^cl)'{G,x+XG,xx) - + 

r 



Px + 2XPxx - 2(G^0 + XG,0jf 
2XG^x (4- - 3*'$' + + - + ) - 4:X^G^xx^' ( - + 



e'^P, 



-(v&' - $') + 'i''^' - 



20'G,^ - + 0'2G,^^ + 20'XG,^x - + *' + $'= 



(6) 
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The continuity equation for the matter fluid is 

P^ + *'(p,„ + P„0 =0. (7) 

This equation can be also derived by combining Eqs. ([3])- ([6]). 

In the following we focus on the weak gravitational background characterized by the conditions j^j ^ 1 and ^ 1. 
Then the dominant contribution in Eq. ^ is of the order of (F/r^)^. In order to make comparisons between each 
term in Eqs. ([3]) and ^ relative to F/r'^, we introduce the following quantities 



eF0 = — p — , £GX = p , eG0 = p , ep 



epx = ^ , = , epm = ^ , (oj 



and 



^ _ P^^</)'r _ XPxx , _ Px4>^'r 

^F(l><l> = p , Apx = — , Apx0 = , 

^ ,(p tX -X 

\ _ XG^xx , _ XG^fjtx , _ G^^^4>'r 
^Gxx = —p^ , '^G<j,x = —p^ , ^G4>4>= — -pr, ■ 

The matter density pm is of the order of (F/r^)$. From the continuity equation ([7]) one has Pm/Pm ^ ^ hi the weak 
gravitational background, so that epm In the following we employ the approximation that all the terms in 

Eq. ([S]) are much smaller than 1. For the consistency with local gravity experiments we require that these quantities 
are at most of the order of $ and ^f. 

From Eqs. ^ and ^ we can express $' and in terms of pm, and (the derivatives of) the field ip. Substituting 
these relations into Eq. ([5]) and neglecting the second-order terms ef relative to Ci, it follows that 

□ * = A*l/Om + M2n0 + , (10) 

where U = —r -\ — , and 

ar^ r dr 

Ml — -rT7 (e^* + 1 - e^*eF0 - e^*eGX - + ep + epx + epm) , (11) 
^ (3-e^^)(eF^ + eGx) 



M3 



1 - 2e^'^ + e^"^ 

4(e2* - 2)ep - (5 - 3e2*)epx + 2(1 - e2*)(e2*ep^ + e^'^eax + ec^) + (3 - e2*)(Ap^^ep^ - 3ep„) 



47*2 

Expanding the gravitational potential $ further and picking up the dominant contributions, we have 



(13) 



Ml - (2 + 6$ - ep^ - eGX - ec^ + ep + e-px + epm) , (14) 

_ ep0 + ecx _ ^0 XG^x ,x 
- 2^-"2P+^^' ^^^^ 

4$^ + Xf<p<p£f<p — 2ep — epx — 3epm 
■ ^1^^ 

Since {|$|,|ez|} < 1, Eq. (Ull) gives mi - 1/(2P). In GR one has F = M^j l/(87rGjv) (Gw is the Newton's 
gravitational constant), so that mi — '^t^Gm- If M2 7^ 0, then the field Laplacian term gives rise to the modification 
to the gravitational constant. This comes from the fact that Eq. ^ contains the matter density pm after replacing 
the term $' by using Eq. ([3]) . Equation ([15]) shows that the gravitational constant is subject to change for the theories 
characterized by 



P^^O, or G,jf^O. 



(17) 
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The right hand side (r.h.s.) of Eq. ^TU\i is of the order of $/r^. Provided that the condition \ii2<l>' /r\ ^ j^aj is satisfied, 
the contribution of the third term on the r.h.s. of Eq. (|10p can be neglected relative to the second term. This amounts 

to the condition {Ief^I, kcxl} > {$^, \>^F<j><j>eF^\, |ep|, \epx\, |epm|}- 

We combine Eq. © with Eqs. (IS])-© to derive the closed-form equation for (p. Picking up the dominant terms, we 
obtain 



where 



□ = Pm + A*5 , (18) 



^ 2F-a ' (^^^ 

/X5 ~ -0'{[2e2*(l - Xgxx) + (e^* - 15)(1 + \GXx)]eGX + 2(26^* - SXq^x + Aq^^ - 2)eG0 

-2(Apxe2* + Apx0 + e^* - 1 - 5Apx)epx - 2ep^}/(2ra) , (20) 



and 



a = (1 + AGxx)(e'* + 3)eGX + 2(1 + \G^x)eG^ - (1 + 2Apx)epx • (21) 
Using the original variables with the approximation e^* ~ 1 in Eqs. ([1^)) and (PU)) . we have 
. „K^,.-'^-/3 + ^-Gx) 

P^r^ + 4X(2G,^x - P,xxWr + [(P^x - G,^,^,^ + 6G,x + ^XG^xxW^ ,oQ^ 

' n > (,20J 

where 

/3 = (Fx + 2XPxx - 2G^^ - 2XG,^x)r - 4(G,x + XG^xxW ■ (24) 
Substituting Eq. (HH) into Eq. HO]), it follows that 

□ * = 47rGoff /Om + + , (25) 



where 



Goff = -^(/il + Ai2M4) 

47r 



1 



SttF 



1 + ( ^ y- I ^ ^ — + 3$ - -(ep0 + ecx + cg^ - ep - epx - epm) 



.(26) 



Equation (P5)l corresponds to the modified Poisson equation. The second term on the r.h.s. of Eq. ([25)) is crucially 
important for estimating the modification of gravity. For the theory in which P — X and G = 0, the second 
term includes the contribution of the order of F'^/(2F). In the presence of the dilatonic coupling of the form 
F = Af^je-^Q^A^Pi, where Q is a constant of the order of 1 and A/pi is the reduced Planck mass, this contribution 
reduces to 2Q^ for |(^/Mpi| ^ 1, so that the gravitational coupling is strongly modified relative to GR. In this case 
the model is in contradiction with local gravity experiments, but the situation is different in the presence of the term 
G{4),X). How this modification works will be the topic of the next section. 



III. VAINSHTEIN MECHANISM 



We study how the Vainshtein mechanism works in the presence of the non-linear field self-interaction G((/), X)[I\(I). 
In doing so, we need to specify the forms of the functions F, P, G at some point. Still, we aim to keep the analysis as 
general as possible, so that it can cover a wide range of scalar-tensor theories. 

First of all, let us consider the form of non-minimal couplings P((/)) with the Ricci scalar R. If we take the power-law 
function P((/)) oc (jf^ the parameter ep^ reduces to tp^ = p(j)'r/(j). Suppose that there is a solution characterized by 
(j) (X . Since ep^ = pq in this case, the approximation jep^j ^ 1 breaks down for p and q of the order of unity. 

For example, let us consider Brans-Dicke theory [s^ characterized by the functions F{(j)) = cj) (i.e. p = 1) and 
P = uj-QTiX/cj), in the presence of the X-dependent term G{X). For those large values of the radius characterized by 
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the condition IPx'^l ^ |4(G,x + XG^xjc )</>'! the term /? in Eq. ([2^ is given by /3 ~ ci;bd?'/0, so that /Z4 is constant 
for the power- law solution (j) cc r''. In fact, as we will see shortly, there is a solution with (j> (x r^^ for constant ^4. In 
the regime l-P.x?"! ^ |4(G,x + XG,xx)4''\ the field equation ([T8)) is approximately given by 

-^(r2 0')-M4Pmr2. (27) 

We introduce the Schwarzschild radius Vg of the source, as 

1 f 

which leads to the relation pm,r^ — M^^drg/dr. Then Eq. (P71) is integrated to give 

</.'(r) ^ . (29) 



For the values of r. for which the function rg is almost constant, one has (f>'{r) cx r^^ (implying pm''^ ^ 0), which 
in turn leads to ep^f, = —1. Hence the validity of the approximation used in Sec. HIl breaks down for the coupling 
F{4>) = (j). This problem can be avoided by rescaling the field (j) with the form of the exponential coupling, i.e. 
(f> — > MpjC^^'^'^/^p', where Q is constant. In this case the rescaled field is related with a dilaton field appearing in 
low-energy effective string theory [s^- As we will see in Sec. lIVBl after this rescaling, the kinetic term in Brans-Dicke 
theory reduces to the standard one P = —g^'^d(j)f^d(j)i,/2 in the limit Q — > 0. We treat the field (p with the exponential 
coupling 

F((/.) = A^e-^Q-^-A^pi (30) 

as a more fundamental one rather than the field (f> with a power-law coupling F{(j)) cx 0^. In this case one has 
eF4> = — 2(50'r/A/pi, so that \eF4>\ can be much smaller than 1 even for the power-law solution like 0'(r) cx r~^. 

In the following we focus on the theories with the coupling ([50)) . In doing so, we shall study two different cases: 
Q 7^ and Q — 0, separately. 



0/0 



When we discuss the case of non-zero values of Q, we are primarily interested in the theories where \Q\ is of the 
order of unity. In this case the term in Eq. (|22p provides an important contribution to the field equation. From 
Eqs. and we find that the qualitative behavior of solutions is different depending on the radius r. The 
behavior of solutions changes at the radius ry characterized by 



\B{rv)rv\ = |4(G,x + XG^xx)irv) <P'irv)\ 



where 



B = P 



X 



2XP 



XX 



2G, 



2XG 



<j,x ■ 



(31) 



(32) 



Here ry is the so-called Vainshtein radius [T^ below which the General Relativistic behavior can be recovered in the 
presence of the term G{(p, X)\I\(I). 

In what follows we focus on the theories described by the action 



S 



-9 



-F{<f>)R + f{^)X - 5(,^)M^-*"X"n0 



(33) 



where F((f)) is given by Eq. ([50]) and M is a constant having a dimension of mass. In this case P{(j), X) = f{(j))X and 
G{(j),X) = f;(0)M^~'*"X", where f{4'),g{4>) are functions of </> and n is a positive integer {n > 1). We assume that 
/((/)) and g{4>) are slowly varying dimensionless functions of the order of unity, such that 



|Mpi/,0//l ^ 1 : \Mpi9,4./9\ < 1 . 
If /(0) and g{(j)) are proportional to F(0) = TV/^je^^Q-A/A/pi^ ^j^g condit ions p4p are satisfied for \Q\ < 1. 



(34) 
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For the action Eq. (HH) yields 

d_ 2 _ r[2QF/Mpi + {/ - 2{n + l)Mi-4"g,^X"}0V - 2n(4n + l)Mi-^"gX"] ^ 
dr^^ ^' ~ 2F[{f - 2{n + l)Mi-4«g^^x"}r - 4n2Mi-4"gX"- V] ^""^ 

J.4,Xr^ + 8nMi-4"g,0X"0'r + {(/^^ - Mi-4"^^^^x")r2 + 2n(4n - l)Mi-4"^x"-i}0'2 
{/ - 2(71 + l)Af _^^«}^ _ 4n2Mi-4"gX"- V ' 



(35) 



The quahtative behavior of the solutions to Eq. ((35|) is different depending on whether r is larger than ry or not. 
Moreover the solution is subject to change for r smaller than r*, where r* is the radius at which the contribution of the 
density dependent term in Eq. psp becomes comparable to the last term in Eq. (|35p around a spherically symmetric 
body. Hence there should be three different regimes: (a) r^ry, (b) r* <C r <C ry, and (c) r <C r*. In the following 
we shall derive the solutions to Eq. ([55)) in each regime. 

1- r'y^rv 

As long as the condition jF^j ^ \<i>' fi\ is satisfied in the regime r ^ ry, the term in Eq. can be estimated 
as ~ Q/{Mp\B). As we will see below, the field behaves as (j>'{r) ~ {QM-p\/ B)(r g / r"^) in the regime r ^ rv for 
B ~ constant. In this case the condition \F^^\ ^ \4>' P\ is in fact satisfied. For the theories given by ([55)) the term 
fJ-iPrn is the dominant contribution to the r.h.s. of Eq. ()18p for r ^ ry. This is known by deriving the solution of 
Eq. ([T5)) without the term /is and by substituting the solution into Eqs. ([^^ and ([251) . Then, in the regime r ^ ry, 
the field equation ([55)) reduces to 

where B — f — 2{n + l)M^~'^"g,^X'^ ~ /. As long as /((/)) is a slowly varying function, B is nearly constant. In 
general, for k-essence theories in which P includes non-linear terms in X |40| . B depends on r. Provided that the 
term of the form P = f{(j))X corresponds to the dominant contribution to B, one can also employ the approximation 
that B is nearly constant. Using the Schwarzschild radius defined in Eq. ([25)) . Eq. ([55)) is integrated to give 



0'(O-^S- (37) 



Here we neglected the solution of the homogeneous differential equation as this is equivalent to the renormalization 
of Tg . Under the approximation that the solution ([57]) is valid at r = ry , it follows that 



'V 



iQMplTg 



{G,x+XG,xx){rv) 



(38) 



For the function G = X/A-fi the Vainshtein radius is known to be ry = \^QMp\rg / {B{ryY M'^)\^/'^ . If the term 
G^x + XG^xx depends on X, we need to use Eq. ([37| again to derive the closed-form expression of ry. 

2. r, < r < ry 

In the regime <C ?' ^ ry the G-dependent terms are the dominant contributions to /3. As long as the solution 
is described by 4>'{r) cx with < p < 1, one can approximate /3 ~ —A{G,x + XG ^xx)<t'' {f)- For the function 
G = g{(j))M^~^"-X"- [n > 1) the solution derived later behaves as (/)'(r) c>c r^^/'^")^ go that the approximation given 
above is justified. In the regime r, <C r <C ry the terms {6G^x + 8XG^xx)4>'^ in the numerator of /is in Eq. ([23)) 
provide the dominant contribution to the field equation, and hence 

d . 3G,x + 4.XG^xx 4n-l 

dr 2{G^x+XG,xx) 2n 

Since Xqxx = XG^xx /G,x = n — 1, the coefficient in front of the term rip' in Eq. ([5^ is constant for the function 
G = g((/))Mi-'*"X"'. Equation ([39]) is integrated to give 

0'(r) = Cr-i/(2n) 
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The coefficient C is approximately known by matching two solutions p7p and at r = ry, which gives C = 
QMpirgry^^"^"V-B(ry). Then Eq. (gO]) reduces to 

,V,.^f^V'""\ (41, 



B{rv)rv \rv 

Compared to the solution ([57)) the field derivative varies more slowly in the regime ^ r ^ ry. This is the region 
in which the Vainshtein mechanism is at work. The solution (j4ip is compatible with the approximations we made to 
find Eq. dai), e.g., \g^^g-^(j>'r\ < 1. 

The solution (|4ip diverges in the limit r — >■ 0. To avoid this divergent behavior we expect that <j)' {r) behaves in a 
different way for the radius smaller than r*. In order to find the radius below which the approximation we used 
breaks down, we compute the variables defined in Eq. ([S]): 

, ^ 2Q^ / r Ml r IQ^%(0)Mg,| ( r^V ( r V-"- 

^ W^WW) w ) w) ' '"-"^^W^WWWJ w) ' ^ ^ 

where ?7/((/)) = Mp\f,^/f and r/g(0) = Mp\g^^/g. As long as F{(f))/Ml, f{(f>), g{(f>) do not change significantly and they 
remain of the order of unity, the quantities given in Eq. (j42p . apart from |eGx|: a-re much smaller than 1 for r ^ ry- 
The variable ecx is proportional to r~^^^'^"'\ which diverges in the limit r — > 0. The validity of the approximation 
I ^GX I ^ 1 breaks down for the radius r < r^, where 



2n 



'^''[AnF{r!)\B{rv)\rl) " ^^^^ 

For \Q\, |i3(ry)|, and F{r^)/Ml of the order of unity one has r^/vg w (rg/ry)*"^^ . This shows that, for n > 1 and 
<C yy, is extremely small even compared to rg. As we will see in Sec. IIVI the typical Vainshtein radius for the 
Sun {rg w 10^ cm) is around ry ~ 10^" cm for the models relevant to dark energy. When n — 1 one has r^ w lO"**" cm, 
which is even smaller than the Planck length. As we will see in Sec. IIVI the typical value of r» is about the radius of 
the Sun (« 10^° cm), so that the solution (41) is trustable for r > r*. 

3. r ^r. 

Let us derive the solution to the field equation in the regime where the density dependent term on the r.h.s. of 
Eq. (|35p becomes important around the spherical symmetric body. For the regularity of solutions the boundary 
conditions should satisfy 0'(O) = and |0"(O)| < oo at the origin. These two conditions lead to (/!)'(r) cx r™ (to > 1), 
as r — > 0. We also impose that the density approaches a constant value pc in the limit r — >■ 0. 

For the theories with n = 1, i.e. G = g{(t>)M^^ , Eq. ([55]) around r = reads 

A(.V) "^^^ ""-^ r (44) 

where fc = f{(f>c) and gc = g{4>c) with (pc being the field value at the origin. On using pm — pc — constant, there is a 
solution characterized by (/''(r) = hr. The coefficient h is known by substituting (t)'{r) = br into Eq. I^^ . This leads 
to the following solution 



(t>{r) 



^9c 



1±W1 --^rcgc 



(45) 



If the condition \Qpcgc\ ^ Mp\M^ f"^ is satisfied, Eq. (|45p is approximately given by 

( \Q\M^I 

V 6Mpig, 



8 



whose existence requires that Q < for (^c > 0. The sign of Eq. is fixed by matching this solution with the one 
in the regime r ^ r». When B{rY) > and B{rY) < the sign of Eq. (pS)) is negative and positive, respectively, for 
g <0. 

While we derived the solution around the center of the star, this is also valid for the star where p„i is 
approximately constant. Matching the two solutions (|1T|) and (|^ at the radius r* for n = 1, it follows that 



6\Q9c\ i V^'^V 



B^{rv)pcrl M 



(47) 



at which the two terms on the r.h.s of Eq. psp are the same order. The matching radius depends on M as well as 
Pc- For the models relevant to the cosmic acceleration today, for the Sun is typically around its radius (provided 
that the density of the Sun is assumed to be nearly constant). 

For n > 1 the field equation satisfying the boundary conditions at r = reduces to 



This is integrated to give 



where we used ^'(0) — 0. In order to match this solution with (|4ip we require that B{rv) > 0. 



4- Corrections to the gravitational potentials 



We estimate the modifications to the Newtonian gravitational potentials in the regime r* ^ r ^ ry. First of all, 
let us see how the Vainshtein mechanism suppresses the additional gravitational coupling appearing in Eq. ([26)1 . Since 
\^Gx \ ^ k-F^I for r^, <^ r <^ rv, the term XG^x/F can be neglected relative to F^^/{2F). Using the solution (j4T|) and 
the definition of ry given in Eq. it follows that \P\ ~ |4(G',x + XG^xxWl - |5(ry)|ry (ry/r)i-i/(2"). Then 

the second term in the square bracket of Eq. (j26p can be estimated as 



2F 



XG 



X 





0''G,x) 








2Q2 




(ry) 






2F/3 




Ml 


B{rv) 





2-l/(2n) 



(50) 



which is much smaller than 1 for r <C J'y- Hence the presence of the term G(0, X) in (3 can lead to the recovery of 
GR within the Vainshtein radius. 

Using the estimation (|^^ and picking up the dominant contributions in Eqs. ([3]) and (jH), it follows that 



2F , 2F 

— $' + — $ - F" 



2F' 

r r 

2F , 2F 2F' 

— + ~0. 



(51) 
(52) 



Substituting the solution (HT|) into Eq. ([?T|) . we obtain 



d 
dr 



(1 - 4n)02 



2nB(rv)4-i/(^") 



9 ^l-l/(2n) 



Integration of this equation leads to 



'_g_ 
2r 



2Q^ 

B{rv) \rv 



2-l/(2ri)' 



(53) 



(54) 



where we neglected the homogeneous solution, as it corresponds to the renormalization of r^. Plugging this solution 
into Eq. ((52|l . we get 



'_a_ 
2r 



An 



2n- 1 B{rv) \rv 



2-l/(2n)' 



(55) 
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Clearly the second terms on the r.h.s. of the square brackets of Eqs. and (IF5|l are much smaller than unity in the 
regime r <C ry, so that the fifth force is suppressed. 
We define the post-Newtonian parameter 7, as 



7 = . 



(56) 



The present tightest experimental bound on 7 is I7 — 1| < 2.3 x 10 ^ ^]. Using the solutions (|54p and ([55]) this 
constraint translates into 



2g2 



1 



2n- 1 \B{rv)\ \rv 



2-l/(2n) 



< 2.3 X 10" 



(57) 



For r much less than ry the bound ([57)) can be satisfied even for \Q\ =0(1). 

In the regime r <C the corrections to the gravitational potentials are even more suppressed than those estimated 
by Eqs. and (155)1 . This comes from the fact that, as r approaches 0, the quantities such as jeF^I and |eGx| 
decrease. 



B. 



Let us study the theories given by the action ^ with F = M^^. Even for Q = there is an additional correction 
term XG^x/F in Eq. ()^^ . so it is not clear whether such a correction is suppressed or not. We discuss two different 
cases: (i) n = 1 and (ii) n > 1. 



1. n = 1 



This case corresponds to the function G = g{(t>)M~^X . The qualitative behavior of solutions changes at the radius 
rv characterized by the condition — gv4>' (rv)\, where /y and gv are the values of / and g ai r = ry, 

respectively. 



For r ^ ry the term ^4,pm dominates over in Eq. (fT5)l . where /i4 ~ 0'r/(2Afp[). It then follows that 



dr 



2 dr 



The solution to this equation can be written in the form 



C 

I2 exp 



9 , 



Pm{r)fdf 



(58) 



(59) 



where C is an integration constant. For the local matter density we assume that the integral J°° Pm{r) f^df is finite. 
Then, for large r, we require ~ br^^^"^, with q > 1 (5 is a constant). In this case Eq. ()59p yields 



(j)'{r) = -^exp 



(60) 



where Ci is another constant which absorbs the exponential term (which is nearly constant for large r). Note that 
this relation is valid even for the weaker bound q > 0. 

Let us consider the regime r < ry. Around the center of the spherical symmetry there should be the change of 
solutions at some radius r*, so we first derive the solution in the regime r* ^ r ^ ry. Since /3 ~ —AAI^^gcf)', 
II4 ~ 5(/)'r/(8Mp[), and /is ~ 3(/)'/(2r) in this region, the field equation is 



5 drg 3 
r I - ' 



.8 dr 2 

Under the condition \drg/dr\ ^ 1 (which corresponds to pmT^ for large r) we have the approximate solution 



(61) 



(62) 
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where C2 is an integration constant. The solution (p^ cannot be trusted up to r — )• 0. 

Finally we study the behavior of the solution in the region r ^ r^. For the regularity at the origin the field 
derivative should take the form (i>'{r) cx r™, where m > 1. Since p„i = 0{r'^) in this regime, the field equation 
psp is approximately given by 



Assuming the solution of the form (j)' = br and substituting it into Eq. , we find the following solution 

^9c 



(63) 



(64) 



Let us match Eq. (|64l) with Eq. ([62|) at the radius r = r*. We caution that, if the density pm differs from p^ 
at r = r*, the correction to the solution (|64|) should be taken into account. Since Ci ^ fc/i^gc) Af3 rl^^ after the 
matching, the solution in the regime rg < r <C is 



*'M = ^"=-(7)"'- 

Finally we match Eq. (|65p with Eq. (|60p at r ~ ry. This leads to the following solution in the regime r 3> ry: 

„ / X 1/2 , 



(65) 



(66) 



The Vainshtein radius ry is defined by the condition l/yryj ~ |4M ^gyip' {rv)\- This gives 



rv 



fvQc 



2/3 



(67) 



For I /I and | g| of the order of unity, Eq. ([57)1 implies that ry is the same order as r*. This means that there is no 
intermediate regime < r < ry characterized by the solution (j65|) . Moreover the matching radius ry cannot be fixed 
completely. 

Let us consider the regime r < ry with the solution (|64)) . Then the second term in the square bracket of Eq. ((26|) 
can be estimated as 



2F F 



(68) 



If the same model is responsible for the late-time cosmic acceleration, the mass AI is related with today's Hubble 
parameter Hq as Af^ w AfpiiJp. Using this relation one has ^ w {t/Hq^)^, which means that the correction is 
significantly suppressed on solar-system scales. The G-dependent terms appearing in Eqs. ([3]) and Q are the orders 
of {F/r^)€Qx and {F/r'^)eQ^. For the radius r < ry we have lec^l ~ (f/HQ^)^, which is even much smaller than 
kcxl- Then the corrections to the gravitational potentials coming from the term G{(l), X)D(f) are strongly suppressed 
on solar-system scales. 

Li the regime r > ry the field derivative (f>'{r) is a decreasing function with respect to r, so the correction ^ becomes 
maximum around r = ry. As long as ry <C Hq^, ^(ry) (ry/_ff(^^)'* is much smaller than 1. Of course, if r,(~ ry) 
is significantly away from the origin, we need to take into account the correction to the solution (j64p coming from the 
change of the matter density. Still the small cubic mass term appearing in Eq. ([M)) would affect the solutions in 
the regime r > ry [as it happens in Eqs. (|65|) and (j66p ]. so that the term ^ should be suppressed as well. 



2. n> 1 



We proceed to the case in which n > 1. Let us study the behavior of solutions around the origin. For the regularity 
we need to assume the form (/!)'(r) cx r™ (to > 1) and pm Pc, as r 0. Then Eq. ([55)) yields 



dr^""^^ 2Af2/<^ 2/(</),) 



(^2^.)^_^^^3^/_^|V^^2^,2_ (gg) 
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For the theories in which / is constant the solution to Eq. ([M)) is given by (j)'{r) oc eP""- /(«^pi)/^2^ ^^^^^ 

is singular 

at r = 0. If / depends on ip, we obtain the following solution 

which is again singular at r — 0. In both cases the solutions cannot satisfy the regularity condition at the origin. The 
theories with n > 1 and Q = are not viable because of the above mentioned property. 

IV. APPLICATION TO CONCRETE MODELS 

In this section we apply our formulas given in Sec. Illll to a number of concrete models. 

A. Extended Galileon 

We first study the theories characterized by 

F{(l)) = M^ie-^'S'^/*^!" , P{X) = eX , G{X) = AM^-^^X" , (71) 

where e = ±1, n is a positive integer {n > 1), and A is a constant of the order of unity (which can be either positive 
or negative). In this case the functions /(</)) and g{(j)) in Eq. (|33p are strictly constant, i.e. f{4>) = e and g{(j)) = A. 

The covariant Galileon modeljwhich recovers the Galilean symmetry dfj_(f> + in the limit of Minkowski 

spacetime, corresponds to n = 1 |28l.[29|. Note that in the DGP model the field self- interaction of the form XM^^XDcf) 
arises from a brane-bending mode. For general n the background expansion of the Universe is the same as that of the 
Dvali- Turner model [42]. When e = —1, Q — 0, and A > there is a de Sitter attractor along which (p ~ constant. 
If this solution is responsible for the cosmic acceleration today, the mass AI is related to the today's Hubble radius 
r, = H^^ w 1028cm via M « (MpV^"r2")i/(i-4") ^^\^. 

Let us consider the case Q 0. Using Eq. ([57]) . the Vainshtein radius ry defined in Eq. (pS)) reads 
. ro3-«,2|,.i/(4n-i) (|Q|Mpir,)(2"-i)/(^"-i) ^ (|g|Mpir,)(2"-i)/(4"-i) 

rv = {2 n|A|)/^ . (72) 

If the mass M has an approximate relation M w (Mpj^^'V^")^/'^^^'*") (as in the case of Q = 0), one has ry w 
(|Q|rg"~^r^")^/^''"^^^. When n = 1 this reduces to ry ~ (IQks''c)^^'^' which recovers the Vainshtein radius ry « 
(rgr2)i/3 in the DGP model for \Q\ = 0(1). For the Sun (r^ w 10^ cm) one has ry w 10^° cm for \Q\ = 0(1). 

From Eqs. (|37p and (|41[) the solutions in the regimes r 3> and r* <C r <C ry are given, respectively, by 
(/)'(r) ~ QMpirg/{er^) and 0'(r) ~ QMpirg / {ery){r /ry)^^/^^^^\ If we consider the case n — 1 and A > with the 
condition |QpcA| > MpiM^, the solution in the regime r < r* is 4>'{r) ~ -[|Q|M^pc/(6AfpiA)]^/2r for e = +1 and 
(j)'{r) ~ [|Q|A/^pc/(6MpiA)]^/2r for e = -1 respectively (where in both cases Q < 0). 

For n — 1, the matching radius r* given in Eq. (|47p can be estimated as 



n ^ f6\QX\ryM^,Y' 

Tg \ PcTy rg M3 j 



(73) 



If the mass M is related to r^ via « -^pir;:^ then the ratio yields r^/rg « [2|QA|(ry/rg)(po/pc)(rc/ry)*]^/^, 
where po ~ 3Mpj/r^ k, lO^^^g/cm'^ is the cosmological density today. For the Sun {pc ~ lO^g/cm^) one has 
r* « lO^rg « 10^*^ cm for \Q\\ — 0(1), which is the same order as the radius of the Sun. The distance r^ given in 
Eq. (|43p is very much smaller than r*, so that the solutions derived in Sec. Illll are trustable. We caution that around 
r = r» there is a correction to the solution (1461) coming from the varying matter density, but still cannot be smaller 
than Tg. 

When n > 1, matching two solutions in the regimes r ^ r, and ^ r <C ry requires that B{ry) > and hence 
e > 0. However the existence of a late-time de Sitter solution (0 = constant) requires e < 0, which means that there 
are no solutions with appropriate boundary conditions around the origin. While the model with positive e may be 
irrelevant to dark energy, it does not possess the discontinuous behavior. 
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If n = 1, the gravitational potentials and ([55]) in the regime ^ r ^ are given by 



$ ~ -2- 

2r 



l-^(- 



r 
rv 



3/2' 



2r 



4Q^ 
e 



3/2' 



(74) 



When \Q\ < 1 the experimental bound (|57p on the post-Newtonian parameter 7 is satisfied for r < 5 x 10^ ry. If 
the relation « Afpir"^ holds, this bounds translates into r < 10^^ cm for the Sun and r < 10^^ cm for the Earth. 
Note that for the radius r ^ the corrections to $ and ^ are much smaller than those given in Eq. ((74)) . Hence the 
model is compatible with the experimental bound on the solar-system scales. 

When Q = we showed in Sec. IIII Bl that the corrections to the gravitational potentials are extremely tiny for 
n = 1, so that the model can pass the solar-system constraints. The models with n > 1 are plagued by the problem 
of the singularity of 0' (r) at the origin. We note that this situation may change in the presence of other non-linear 
field corrections appearing as the forms of £4 and £5 in Galileon gravity [28L [2g | . 



B. Brans-Dicke theories with a universal dilatonic coupling 



The Brans-Dicke (BD) theory fsS*] is characterized by the action 

^BD 



S 



'9 



2X 



(Vx)^ 



(75) 



where x is the scalar field coupled to i?, and ojbd is the BD parameter. Introducing the field in the form x ~ 



M. 



pi^ 



-2Q0/A/pi^ the action ([73) can be written as [14] 



where 



1 



2(3 + 2a;BD) 



(76) 



(77) 



If we define the field ip as ip = 2(50, the square bracket in Eq. (|76|) is expressed as £ = M^^e v/Mpi [/j_[jgQ(v^)2]y'2-|^ 
• • • . This means that dilaton gravity [39| corresponds to ojbd = ~1- Motivated by dilaton gravity, we shall consider 
the theories in which the field (p has a universal coupling F{(l)) — with XD(f> as well: 



F{(b)R + {l-6Q^) 



Ml 



X - 



XDcj) 



(78) 



where A is a constant of the order of unity. The last term appears as the a' correction in low-energy effective string 
theory [ij. We assume that the coupling Q is of the order of unity with 7^ 1/6. The action ((78)) corresponds to 
the theories with /((/)) = (1 - 6Q'^)F{(t))/Ml, g{(j)) = \F{(j))/Mly, and n = 1 in Eq. ((33)). We consider the case in 
which the field satisfies the boundary condition 1 0(0)1 ^ Mpi at the origin. Since the field derivative (/)'(r) is small, 
the condition \<t>{r)\ ^ Mpi is satisfied for r > 0. From Eq. ((77)) dilaton gravity (wbd — —1) corresponds to ~ 1/2. 
From Eq. ((55)1 the Vainshtein radius is given by 



TV 



|4QA| Mpirg N 



1/3 



(79) 



,(1-6Q2)2 M3 ; 

where we used the approximations F « and i? 1 — QQ^ . If the model (1751) is responsible for the late-time cosmic 
acceleration, one can show that there is a de Sitter solution characterized by = constant with k, Mpir^ 



For 



r2)l/3. 



\Q\ and |A| of the order of unity the Vainshtein radius ((7^ is estimated as ry sa (r 

Under the condition \QpcX\ > MpiA'P{l - eQ^)^, the solut ion in the regime 
-[|Q|AfVc/(6MpiA)]i/2r for < 1/6 and 0'(r) ~ [\Q\M^pc/i6Mpi\)]^/^r for > i/g. 
r* given in Eq. ((47| can be estimated as 



2\QX\ rv po ( r, 



(1 



^^)^ Tg Pc \rv 



-,1/3 



is given by (f>'{r) ~ 
The matching radius 



(80) 



13 



where we used the relations « Mpir~^ and po ~ 3Mp[/r^. The radius is the same order as Eq. (|73p for 
\QX\ = 0(1), so that the solutions to the field equation derived in Sec. Illll are again trustable. 
Equations ([54|) and ((55|) show that, in the regime ^ r ^ ry, the gravitational potentials can be estimated as 



2r 



1 - 



1-6Q2 



3/2' 



2r 



4Q2 



i-6g2 



3/2' 



(81) 



The experimental bound (|57p translates into 



3/2 

— I < 2.3 X 10" 



|1-6Q^ 
2Q2 



(82) 



In dilaton gravity (Q^ = 1/2), for example, this constraint is satisfied for r < 10 ^ry ~ 10"'^^ cm for the Sun. This 
upper bound is much larger than the solar-system scales. 



V. CONCLUSIONS 



In this paper we have studied the Vainshtein mechanism in second-order scalar-tensor theories given by the action 
([1]) . In a spherically symmetric background the full equations of motion were derived in the presence of a barotropic 
perfect fiuid. Introducing the small parameters defined in Eq. ([5]) and picking up the dominant contributions in 
the weak gravitational background, we obtained the closed-form equations for the field cf) and for the gravitational 
potential ^. The approximation employed in Sec. |n]is valid under the conditions je^l ^ 1, which is required for the 
consistency with solar-system experiments. 

The general theories in which the Vainshtein mechanism can be at work are given by the action (|33p with the 
nonminimal coupling F[(j)) = Mpje~2'30/A^pi , This action covers a wide range of modified gravitational theories such 
as the (extended) Galileon, dilaton gravity, and Brans-Dicke theories with the nonlinear field interaction. In such 
theories we derived the general formula ([55)) for the Vainshtein radius ry- 

For Q ^ the solution to the field equation in the regime r ^ ry is given by Eq. ([57]) . which leads to the large 
modification to the gravitational potentials. In the regime ^ r ^ ry the solution changes to Eq. (|41l) . so that the 
modification of gravity is suppressed even for \Q\ — 0(1). In this regime we derived the analytic solutions for the 
gravitational potentials and showed that the experimental bound on the post-Newtonian parameter 7 can be satisfied 
under the condition ([57]). In the regime r ^ the solution to the field equation for Q 7^ is given by (f>'{r) oc r, 
which satisfies the regularity condition 4>'{0) = at the center of the spherical symmetry. 

When Q = and n = 1, if the same model is responsible for the late-time cosmic acceleration, the corrections 
to the gravitational potentials are extremely tiny on solar-system scales. For Q — and n > 1 we showed that the 
solution around the origin is not regular and hence this theory cannot be regarded as a viable one. 

We applied our general results to concrete theories such as extended Galileon and Brans-Dicke theory with a 
dilatonic coupling. For the theories with Q ^ and n = 1 there is a correction of the order of (r/ry)^/^ relative to 
the Newtonian gravitational potentials in the regime r ^ ry-, such that the local gravity constraints can be 

satisfied on solar-system scales. Note that for Q = and n — 1 the corrections to the gravitational potentials are 
even much smaller. If the extended Galileon theory with Q 7^ and n > 1 is responsible for the cosmic acceleration 
today, we found that there is a problem of the matching at r = . 

It will be of interest to see how the Vainshtein mechanism works in the Horndeski's most general scalar-tensor 
theories having the term G5 as well as the X-dependence in G4. The construction of viable dark energy models 
satisfying recent experimental and observational bounds will be also interesting. In particular the constraint coming 
from the variation of the Newton constant {\G/G\ < 0.027Jo) can provide tight bounds on such models j45j. We leave 
these issues for future work. 
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